PROCEEDINGS 

OF SCIENCE 



DESY 10-019 
SFB/CPP-10-21 



o 

(N 
IX, 



Status and prospects for the calculation of 
hadron structure from lattice QCD 



in: 



Dru B. Renner 

NIC, DESY, Platanenallee 6, D-15738 Zeuthen, Germany 
E-mail: |dru . renner@desy . de 



> 

in 

(N 

o 

o 
o 



Lattice QCD calculations of hadron structure are a valuable complement to many experimental 
programs as well as an indispensable tool to understand the dynamics of QCD. I present a focused 
review of a few representative topics chosen to illustrate both the challenges and advances of our 
community: the momentum fraction, axial charge and charge radius of the nucleon. I will discuss 
the current status of these calculations and speculate on the prospects for accurate calculations of 
hadron structure from lattice QCD. 



X 



The XXVII International Symposium on Lattice Field Theory 

July 26-31, 2009 

Peking University, Beijing, China 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Status and prospects for the calculation ofhadron structure from lattice QCD 



Dm B. Renner 



1. Introduction 

Our most stringent constraints on the structure of any hadron follow from the underlying sym- 
metries of QCD. Translational invariance dictates that the momentum of a hadron is the sum of 
the momenta of all its constituents, giving rise to a powerful sum rule for the nucleon. Rotational 
invariance demands that the spins of hadrons have the values n/2 for an integer n, hence the nu- 
cleon spin is precisely 1/2 and not 0.5 with some experimental error. Similarly, we know that the 
electric charge of the nucleon is exactly 1 and the net strangeness is 0. These statements are so 
commonplace that they may even appear trivial, but each of these exact results serves as an entry 
to different aspects of the dynamics of QCD. 

The momentum of the nucleon is built up from that of the quarks and gluons. While trans- 
lational symmetry constrains the entire contribution, the individual contributions of each quark 
flavor and the gluon depend on the details of QCD. Similarly, the spin of the nucleon arises from 
the quark spin, quark orbital motion and the gluon angular momentum, each of which are individ- 
ually unconstrained by symmetries but must sum to 1 /2. The charge of the nucleon arises from the 
quark charges that yield a total charge of 1 , but the distribution of this charge in the spatial degrees 
of freedom probes the nonperturbative structure of the nucleon. Finally, the strange quarks in the 
nucleon occur in precisely matched pairs of quarks and anti-quarks. Despite the net absence of 
strangeness, the consequences of this hidden flavor are felt in many observables. In the following, 
I will elaborate on a few of these points and discuss what we currently know from lattice calcu- 
lations, phenomenology and experimental measurements. Along the way, I'll offer my opinion 
regarding where our calculations may have to go to provide accurate results for hadron structure. 

The past several years have seen extensive reviews of the calculations presented each year 
at the annual Lattice conference. You can find the few most recent reviews in ||, ^j], and a 
very recent and exhaustive collection of results can be found in [ffl]. Rather than duplicating these 
efforts, I will instead present some of the key examples mentioned above. 1 This will unfortunately 
prohibit me from discussing all the hadron structure efforts presented at Lattice 2009, but I hope 
these proceedings will still provide a useful overview, nonetheless. 

2. Nucleon Momentum 

The momentum of a nucleon with energy E and mass m is precisely constrained by Lorentz 
symmetry as p 2 = E 2 — m 2 . This is such a common statement in particle physics that it seems nearly 
trivial to even mention it. However, here we want to question how this momentum arises from the 
underlying QCD degrees of freedom. In other words, what is the distribution of this momentum 
among the nucleon's constituents? This question concerns the details of the dynamics of QCD and 
presents a challenge to lattice calculations of nucleon structure. 

2.1 Parton Distribution Functions 

The proper field theoretic response to the question above requires constructing the momen- 
tum distributions of the nucleon's constituents, the parton distribution functions (PDFs). I'll pro- 

1 At the lattice conference I discussed strangeness in the nucleon. A recent review [Q] covers much of this, so I will 
not include these results in this review. 
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Figure 1: Momentum distribution of quarks and gluons in the nucleon at jj. = 2 GeV. For each parton 
xf(x,)J.) is plotted where f(x,fl) is the corresponding PDF. At large x, the up and down quarks are the largest 
components of the nucleon momentum. At low x, the gluon dominates and the anti-quark distributions grow. 
The curves were generated using the LHAPDF library [g] and the MSTW2008 NNLO [0] dataset at a 
renormalization scale of = 2 GeV. 

vide a proper definition of the PDFs shortly, but for the moment we want to focus on a more 
intuitive understanding of the quark and gluon distributions, q(x) and g(x). In the parton model 
q(x) or g(x) would be the probability to find a quark, of some flavor, or a gluon with momentum 
^parton = ^nucleon- This interpretation is retained in QCD, however, q{x,jx) dxi&g(x,fx) now carry a 
renormalization scale pL that loosely gives the energy resolution at which the distribution is probed. 
While this muddies the picture a bit, the PDFs, nonetheless, are well-defined universal properties 
of the nucleon that are probed in many different experiments. 

2.1.1 Phenomenological Results for PDFs 

There is now a well-established industry dedicated to extracting the PDFs from global anal- 
yses. The analyses differ in many details, including the order of the perturbative expansion used, 
the treatment of the quark masses and the handling of the statistical and systematic errors from the 
many experimental inputs. These variations have an impact on the precision ultimately obtained, 
however, these details do not impact the discussion at hand. Hence, as merely one example among 
several, I show the PDFs from the MSTW collaboration in Fig. [jj. 

Field theory effects ultimately handicap any simple interpretation of the PDFs, however, we 
can still see that the PDFs at a low scale, /i = 2 GeV in Fig. [T], retain features that one might expect 
for the nucleon. The up and down quark distributions are the largest component for high momentum 
(x < 1). These distributions in fact have peaks at reasonable values of x, as one might expect for a 




3 



Status and prospects for the calculation ofhadron structure from lattice QCD 



Dm B. Renner 




Figure 2: Momentum distribution of quarks and gluons in the nucleon at ji = 100 GeV. The details are the 
same as in Fig. |l] but now fx = 100 GeV. At high scales, the PDFs for the quarks and gluons mix and the 
simple non-relativistic nucleon structure, still visible at jx — 2 GeV, is suppressed at pi = 100 GeV. 



hadron dominantly composed of two up quarks and one down quark. However, the most prominent 
feature in the plot for low x is clearly the gluon distribution. This is a clear indication of QCD 
physics at work. This has a counterpart in the simple observation that the quark masses directly 
contribute only about 1% of the nucleon's mass. Beyond the dynamics of QCD, the presence of 
the anti-quarks is a striking field theory effect. This statement may seem mundane, but we must 
remember the role played by the anti-quarks in nucleon structure to appreciate why we go through 
all the difficulty to calculate hadronic matrix elements from fully dynamical lattice QCD rather 
than being content with models or quenched calculations. 

As mentioned, field theory complicates the interpretation of the PDFs. To illustrate this, I 
show the same distributions again in Fig. |2]but now evolved to a scale of /I = 100 GeV. Any hint 
of nucleon physics is now well hidden and all the constituents of the nucleon appear to play nearly 
equally important roles. We will return to this issue again when discussing the evolution of the 
momentum fraction. 

2.1.2 Operator Definition of PDFs 

The parton distribution functions can be defined as nucleon matrix elements of quark and gluon 
fields separated by light-like distances. As an example, we record here the operator definition of 
the unpolarized quark distribution q(x,Li 2 ) 

q(x,Li 2 ) = IJ^ e u ^ n (p,s\ q(-X/2 n) jt W n {-X/2 n,X/2 n) q(X/2 n) \p,s)\p . (2.1) 
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For each flavor there are two other twist-two quark distributions, the helicity and transversity dis- 
tributions, denoted as Aq(x,/jL 2 ) and Sq(x,n 2 ) respectively. The factor W„(— A/2 n,X/2 n) is the 
Wilson line extending along the arbitrary light-cone direction n from A/2 n to —A /2n, 

( i Xl1 \ 
W„(-A/2 n, A/2 n) = &exp ig I daA(an) - n) . 

V ./-A/2 / 

The scale, fi, and scheme dependence of q(x,jj. 2 ) comes from the renormalization of the operator 



in Eq. 2. 1 



The expression in Eq. ^T] highlights the difficulty of lattice calculations of PDFs. Lattice 
QCD calculations are performed in Euclidean space, however, the PDFs involve quark and gluon 
fields separated along the light-cone. These inherently Minkowski-space observables are difficult 
to construct explicitly in Euclidean space. But, as we will see in the next section, the moments in x 
of these distributions can be calculated directly in Euclidean space. 

2.2 Moments of Parton Distributions 

As just discussed, the PDFs as a function of x are essentially Minkowski-space objects. But 
the moments in x are related to local operators that can be calculated in Euclidean space on the 
lattice. 

2.2.1 Mellin Transform: from x to n 

The moments in x are defined as follows, 

= J dx^qix,^ 2 ) = J dxx" {q(x,ll 2 ) - (-\) n q(x,ll 2 )} . 

The sign in the above equation is determined by the identification of q{x) for negative x with the 
anti-quark distribution as q(x, /I 2 ) = —q(—x,/J, 2 ). 2 These moments can be related to forward matrix 
elements of twist-two operators. First we present the complete result but then sketch the argument 
in the following. 

(p,s\ q(0) iD* • • • iD^qiO) \p,s)\^ = 2<x">^ 2 p^ . ..p^ (2.2) 

The brackets in r^ 1 ^ denote symmetrization of the indices of the tensor T and subtraction of 
the traces. The precise meaning of this operation will be clarified shortly. 



The derivation of Eq. ^2| is almost elementary, but it is so essential to the method underlying 
the lattice calculations that we sketch the arguments using the unpolarized PDFs as an example. The 
first step is to introduce light-cone coordinates and gauge. This is a useful first step in understanding 
several aspects of the PDFs. In defining the light-cone coordinates, you introduce two light-cone 
directions rv\ = (l,0,0,±l)/v2- The two new light-cone coordinates are given by v • n± = v T 
for a four-vector v. (Similarly, ^ , = y • n± = y^.) Light-cone gauge is the choice of A(x) ■ n = 0, 



2 This is a frequent source of confusion when comparing to phenomenological determinations of the PDFs. Also note 
that the sign is different for the transversity distribution but the same for the helicity distribution: 8q{x, jl 2 ) = 8q(—x, jJ. 2 ) 



and Aq(x,[l z ) = -Aq(-x,fl 



2\ 
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which sets the Wilson line W n (— A/2 n, A/2 n) to 1. Then choosing n = «_ in Eq. 2.1, imposing 
the light-cone gauge for ra_ and relabeling A as y~ gives 

q(x,pi 2 ) = l -J^- e ^ + y-{p,s\ q(-y-/2) y+ q(y~ /2) \p,s)\^ 2 . (2.3) 

The next step in relating the moments to local operators is to use the known limited support of 
q(x) to the region — 1 < x < 1 to expand the integration range used to define the moments, 

/l roo 
dxx n q(x,pl 2 ) = J dxx"q(x,pl 2 ) . 

The remaining steps are relatively elementary. The key sequence follows from combining the above 
with Eq. O 

/CO p CO 

dxx n e ixp+r = {ip + Y n / dx(d+) n e ixp+r = 
-CO ,J — CO 

= {-l) n (ip + )- n r dxe ixp+y ~ (d + ) n = {p + )- n 8{p + y-){id + ) n 

J— oo 

where d + = d/dy~. The above manipulations must be understood as acting under the / dy~ in 
Eq. 23. The final result is 

<*"W =2- Vr ( " +1) (/M| ?(0) 7 + (id + ) n q(0) \p, S }\^. 

By inspection, this expression can be seen as the light-cone coordinate, light-cone gauge form of 
the following 

(p,s\ q(0) 4 {in-D)"q{0) \p,s)\^ = 2(x")^ 2 (n ■ p) n+l . 

This expression relates (x n ) q „2 to diagonal matrix elements of local operators. The more familiar 
form follows from writing the above as 

n^...n^{p,s\ q{0) f iD^ ■ ■ -iD^q(0) \p,s)\^ 2 = 2(x n ) q ^ 2 . . .n^p* . . .p* . 



This form makes it clear precisely what the symmetrization and trace removal in Eq. [Lj means. 

2.2.2 Phenomenology of (x) u _d 

The global analysis illustrated in Figs. [I] and ^ can also be used to examine the moments of the 
PDFs. This time taking as an example the results from the CTEQ collaboration, I plot the results 
for (x) q and (x) g in Fig. ^} These results follow from numerically integrating curves similar to 
those in Figs, [j] and ||. Furthermore, the state-of-the-art results for the particular quantity of interest 
to us, {x) u -d, are collected in Tab. [j]. The phenomenological analyses often present results in terms 
of the so-called valence distribution q v (x) = q(x) —q(x) and the anti-quark distribution q(x), but 
Tab. [I] gives the correct combination that is comparable to lattice results. 

What is particularly interesting about the momentum fraction (x) is that it obeys a sum rule, 



3 The error given in is asymmetric, but here we take only the upper error to compare to the lattice results that 
are all higher than the phenomenological values. 
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Figure 3: Running of the momentum fraction of the quarks and gluons in the nucleon. As // increases, the 
gluon momentum fraction (x) g increases towards the asymptotic value of 4/7 and (x) q , for each of the quark 
flavors, approaches the common value of 3/28. (These values hold for the 4 flavor theory.) Despite the 
unique limiting values for x, the non-perturbative input at low jj. clearly dictates the values of the momentum 
fraction for each parton over many orders of magnitude. The results come from the LHAPDF library [|] 
using the CTEQ6.6C [f| dataset. 



where Y* q 1S ine sum over all relevant flavors. This sum rule imposes constraints on the scale 
evolution of (jc}„2 such that (x)„2 asymptotically approaches a perturbatively calculable limit for 
large jx. 

lim (x) 2 

II — ^oo 1 ~ 



16 + 3N F 

lim (x) 



pj^v'/ftii i6 + 3N F 

For Afp = 4, we have lim (x) q = 3/28 10% and lim (x) g = 4/7 60%. However, we can clearly 
see in Fig. || that the asymptotic results for (x) bear little resemblance to the results at any reason- 
able value of pi. In fact, the hierarchy shown in Fig. || is quite clear. The gluons carry about 40% of 
the nucleon momentum, the up and down quarks carry about 35% and 20%, and the strange takes 
most of the remaining 5%. It is precisely this pattern that we eventually hope to understand from 
lattice calculations. 

2.3 Lattice Calculation of (x) u a 

The lowest non-trivial moments of the unpolarized quark PDFs are {x) q . 4 These are the quark 



contributions that enter the momentum sum rule in Eq. [2.4| . Here I use the particular combination 



4 The lowest moments (1}„ are known in terms of the quark valence structure of the nucleon. 
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l x ] Uv -d v 


(x)u-d 


ABMK 


0.1790 ±0.0023 


0.1646 ±0.0027 


BBG 


0.1747 ±0.0039 


0.1603 ±0.0041 


JR 


0.1640 ±0.0060 


0.1496 ±0.0062 


MSTW 


0.1645 ±0.0046 3 


0.1501 ±0.0048 


AMP06 


0.1 820 ±0.0056 


0.1676 ±0.0058 


BBG 


0.1754 ±0.0041 


0.1610 ±0.0043 



Table 1: Phenomenological values for at p = 2 GeV. For each calculation we give the moment of the 

non-singlet valence distribution, denoted by [x] u _ d . (The square brackets denote simply the integral over 
the region < x < 1 as opposed to the region — 1 < x < 1 used in the angular brackets.) These values were 
coUected in [|k1 The original references are [0 (ABMK), [[TJ ||| (BBG), [[i]| (JR), (MSTW), [|B]| 
(AMP06) and [llj [l|] (BBG N 3 LO). This is combined with the result [x]- _- { = -0.0072 ± 0.0007 from 
[p~0|] to produce a result for (x) u -d, which can be compared to lattice calculations. 



{x) u -d as a benchmark observable to evaluate the lattice calculation of moments of PDFs. 

(p, S \qy^iD v h 3 q\p, S ) = 2(x) u _ d ^p^ 
p. 2 

Here q = (u,d) is the doublet of the light quarks. The flavor combination u — d eliminates discon- 
nected diagrams, which would otherwise require a substantial computational investment. Addition- 
ally, this combination also eliminates any mixing with gluonic operators, thus greatly simplifying 
the renormalization of this quantity. 5 Without the mixing, this observable then only requires a 
multiplicative renormalization. Finally, you can calculate the bare matrix elements using nucleons 
with p = 0. Hence, (x) is basically the most accurate scale-dependent observable in nucleon 
structure that we calculate on the lattice. 

Figure |] illustrates all the dynamical lattice QCD results for (x) u -d for pion masses less than 
700 MeV. The values for (x} u -d, m n and a come from a variety of already published sources and 
numerous private communications. (The references are provided in the caption to Fig. The 
most striking feature that you observe in Fig. |] is that, despite the community's efforts to calculate 
with many actions, several lattice spacings and volumes and a broad range of pion masses now 
approaching 250 MeV, the calculations still overestimate the experimental measurement by at least 
30% and maybe as much as a factor of 2. The spread amongst the groups obviously suggests some 
systematic variations, and I will examine two possible explanations shortly, but it is important to 
note that the individual calculations are visibly much more consistent with themselves than with 
the other calculations. Additionally, notice that some groups are beginning to perform calculations 
"at the physical point," as illustrated by the QCDSF calculation in Fig. f|. This phrase quite often 
means simply m n < 200 MeV. Nonetheless, the next generation of lattice calculations will likely 
shed some much needed light on the chiral behavior. However, the discrepancy amongst the groups 
at pion masses where we should be able to reliably calculate today is an issue that needs to be 



5 Strictly speaking, the disconnected diagrams and the mixing with gluons only vanish for lattice actions with an 
exact flavor symmetry. 
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Figure 4: World's dynamical lattice QCD results for The lattice calculations all overestimate 

the phenomenologically determined value. They also show a fair amount of scatter amongst themselves. 
The possible role of finite size and renormalization effects are discussed in the text. The lattice results 
from @ 



17 



HJ (RBC N F = 2 + 1), [|19j [20p (RBC N F = 2), [|2j (LHPC), [22] |2j (ETMC) and 
24p (QCDSF). The experimental result is generated using the LHAPDF library [§] and the CTEQ6.6C 
dataset M. 



addressed. Without resolving these discrepancies, it will be hard to confidently establish physical 
results even with calculations approaching the physical pion mass. 

Before discussing finite size and renormalization effects, I want to make a quick comment on 
the lattice spacing dependence of [x) u -d- This is currently poorly studied. Of the five calculations 
shown in Fig. |], only QCDSF has calculated beyond a single lattice spacing. However, even in that 
case, the range in a that is used to establish scaling is not large. Their results are an encouraging 
hint that lattice artifacts are not a substantial part of the discrepancy in Fig. f|, but there is nothing 
universal about such effects and all the groups must make a stronger effort to calculate at multiple 
lattice spacings. 

A persistent concern in nucleon structure calculations is the role of finite size effects. In 
fact, the results at Lattice 2009 have added much to this issue even if they haven't resolved it. 
In Fig. H|, I examine several finite size studies by various collaborations. Excluding the lightest 
calculations at m n = 260 MeV, one observes no statistically significant finite size effects for any 
of the remaining calculations. These results are consistent with the common rule-of-thumb that 
m n L w 4 is sufficient. 6 However, the recent results of QCDSF at m n = 260 MeV potentially stand 
in contrast to the finite size dependence observed at higher pion masses. This calculation suggests 
that m n L = 4 is at best just barely sufficient to capture the large volume limit of (x) u -d- It is 

6 My use of m n L to gauge finite size effects is, of course, not strictly correct. I am loosely assuming a discussion in 
or near the chiral limit in which 1 jm% will be the dominant length scale. 
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Figure 5: Finite size studies for (x) H _ f /. The dependence of (x)„_j on m K L is shown for several calculations. 
Notice the essentially fiat behavior for m n L > 4 for all but the lightest pion mass. The QCDSF calculation at 
the lightest pion mass may suggest the emergence of a more substantial finite size effect as m K is decreased. 
The curve is a very simple scaling to the m K = 260 MeV results and is meant solely to guide the eye. The 
results are from the same references given in Fig. fjj. 



possible that this observation would vanish with higher statistics. In fact, the three calculations at 
m K = 260 MeV do essentially agree statistically, but the trend in the results is suggestive, especially 
in comparison to all the other results at larger m n . Of course, it is also possible that (x) u -d would 
drop even further with still larger L. 

In a review of this nature it is difficult to perform a detailed infinite volume limit of all the 
results presented at the conference. However, we can illustrate the impact of finite size effects 
by simply making the crude restriction to m n L > 4. This is illustrated in Fig. |6[ The picture 
is certainly clearer and one may even be left with the impression of some downward curvature for 
some of the lattice calculations, but we must guard against wishful thinking. Excluding the QCDSF 
results, each of the groups is statistically consistent with a linear dependence on nv\ with # 2 /dof < 
0.75. The one exception is QCDSF. A linear fit gives # 2 /dof = 2.1. This is not overwhelming 
evidence of non-linearity and all the apparent curvature comes solely from the lightest point at 
m K = 260 MeV. This point is 2.8a less than the second lightest QCDSF point. Increased statistics 
in these calculations and confirmation of this new finite size behavior by other groups will be 
necessary to understand what is happening at m % = 260 MeV. 

Given the results in Figs. ||and ^, there is no concrete conclusion that we can draw yet regard- 
ing finite size effects, but I hope that Fig. Q will stand as a warning that the finite size effects can, 
and likely do, change substantially as we decrease the pion mass and that collaborations pushing 
towards the physical pion mass must keep these effects in mind. At this point one may argue that 
we should investigate the finite size effects as indicated from chiral perturbation theory. It is my 
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Figure 6: Large volume results for (x) u _ c j. The results are from the same references given in Fig. [|, but only 
those calculations with m % L > 4 are shown. The results from each calculation are consistent with a smooth, 
nearly flat m K dependence. However there is a systematic variation between the calculations. 



personal opinion that, given an absence of any real evidence for curvature in {x) u -d an d given that 
such non-analytic behavior is the hallmark of chiral perturbation theory, it is theoretically question- 
able to force a fit to chiral perturbation theory. While this is my opinion, it was clear at the lattice 
conference that many presenters also simply chose to not show chiral fits. Of course, calculations at 
still lighter pion masses may show the expected chiral behavior and systematic errors at the current 
pion masses may obscure this behavior, but we should keep in mind that in the end calculations 
down to the physical point may be necessary to convincingly establish results for (x) u -d. 

One obvious feature of Fig. |6| is that the results from each collaboration have a flat behavior 
but there are clear shifts between the groups. One natural concern in this respect is the renormal- 
ization of the operator used to determine (x)„-d- This is multiplicative and, importantly, quark 
mass independent. However it does depend on the lattice action and hence will vary between the 
different calculations. By taking the ratio of (x) u -d to the value {x) r ^_ d at some canonical refer- 
ence mass, which I arbitrarily choose to be m n = 500 MeV, we can eliminate the renormalization 
factor consistently for each calculation. The result of this is shown in Fig. ^, and we find that the 
various calculations appear to collapse onto a universal curve. As with the finite size effects, we 
can not make a specific conclusion, but this result is suggestive of potential problems in the renor- 
malization of (x) u -d- There are other possible explanations of the systematic shifts between the 
groups, such as the lattice spacing effects mentioned earlier or systematics due to plateaus that are 



too short in the bare matrix element calculations [19, 25]. Independent of the ultimate explanation 
for these problems, I hope that the difference between Fig. |6| and Fig. [7| encourages us to examine 
the systematics of our calculations. 
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Figure 7: Renormalization free ratio: (x) u -dl '{ x Yu-d- The calculations with m K L > 4 from Fig. ^|are used to 
form the ratio {x) u -d I { x Yu-d> where the reference scale is chosen as ;«J f = 500 MeV. This ratio eliminates 
the renormalization factors from each calculation. The calculations are now all consistent with each other. 
This demonstrates that a systematic error, possibly in the renormalization, could possibly account for the 
discrepancy among the collaborations. 



3. Nucleon Spin 

As discussed in the introduction, the nucleon spin is exactly 1/2 due to rotational symmetry, 
and similar to the momentum fraction, it obeys a sum rule [|2£ 



z z q q 

which relates the nucleon spin to the contributions from quark helicity (l) q and orbital angular 
momentum L q and a net contribution from the gluons J g . The asymptotic evolution of the total 
quark contribution J q = | ( 1 ) Aq + L q and the total gluon contribution J g is given by 

1- j - 1 3Nf 

Pl J W 2 2 16 + 3Ar F 

hm J„ .,2 = . 

^ M l6 + 3Af F 

For N ¥ = 4, we find lim J q = 2/7 w 0.60 • 1 /2 and lim J g = 3/14 w 0.40 • 1 /2. Again we find the 
gluons playing a substantial role, now in the nucleon spin and earlier in the nucleon momentum. 
This is a particularly surprising conclusion in light of the naive quark model result £ ? (I)a^ = 1 
and J g = 0. However, similar to the asymptotic results for the momentum fraction, we should be 
suspicious that non-perturbative QCD dynamics can alter the picture at low scales. 
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Unfortunately, much less is known experimentally regarding the decomposition of the nucleon 
spin. The one quantity that is known well is the axial coupling, gA = (1)a« — that is mea- 

sured accurately in neutron beta decay. A recent review [^7|] gives gA = 1.2750(9) and the most 



recent PDG [ )28| ] world average is gA = 1.2694(28). The remaining quark contributions are es- 



sentially the only other known pieces. The experimental result from HERMES in 2007 y29| ] gives 
(1)a« = 0.842(12), (\) Ad = -0.427(12) and (l) Al = -0.085(17). QCD sum rules [g] or model 
estimates [ |3l| ] indicate that J g 0.25 at low scales. Taking the experimental results for (I)a^ and 
this estimate for J g , we can estimate that 50% of the nucleon spin is given by the gluons with the 
remaining divided into 33% from quark spin and 18% from quark orbital motion. This picture is 
much less certain than the momentum sum rule, but again, it provides a challenge to the lattice 
QCD effort. 

3.1 Lattice Calculation of gA 

The nucleon axial charge can be defined by 

or gA = (l)Au-Arf- The moments (I)a^ are the lowest moments of the polarized PDFs. For the 
unpolarized distributions, the moments (l) q are fixed by the valence structure of the nucleon, but 
this is not the case for the polarized distribution. This difference actually makes the moments (1)a? 
even simpler than the momentum fractions, (x) q , discussed in the previous section. In fact, using 
a lattice action with chiral symmetry would eliminate the need to renormalize the lattice operator 
used to determine (I)a^ altogether. 7 However, many lattice actions in use today do require this 
renormalization, but it is not scale dependent and is in many ways simpler than the renormaliza- 
tion required for (x) q . Additionally, since the power of x in (1)a<j is one lower than in (x) q , the 
lattice operator does not contain a derivative and is a simple quark bilinear. Consequently, the 
bare matrix elements can be calculated more accurately for (I)a^ than for (x) q . As in the case 
of (x} u _d, we focus on the u — d combination again to eliminate the computationally demanding 
disconnected diagrams. Furthermore, this combination gives the axial coupling gA = (I)ah-aj that 
is very accurately measured in neutron beta decay. 

In Fig. H| I collect all the full QCD lattice results for g A , again with m n < 700 MeV. Similar 
to (x) Ll ^ ( j, gA serves as a benchmark observable for lattice calculations of nucleon structure. As 
in Fig. |4] for (x) l( _d, we find that the lattice results show a degree of scatter and all consistently 
underestimate the experimental measurement. However, the discrepancy, between 10% and 20%, 
is more mild and the scatter in the lattice results is less severe. In fact, as I try to argue next, it 
appears that this scatter may be almost entirely accounted for by finite size effects, at least to the 
current level of statistical precision. 

In Fig. ^ I examine several finite size studies that are available for gA- The results for the 
heaviest values of the pion mass show a strong finite size effect. This has led to the current view 
that m n L > 6 may in fact be needed to reliably determine gA- However, we find that, as quantified 
by m n L, the finite size effects are diminishing as m n is lowered. This can be seen quite clearly in 



7 This is only true for the non-singlet moments, such as (1)a«-aj- The axial anomaly generates mixing for the 
singlet moments. 
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Figure 8: World's dynamical results for g A . The lattice results are from |J32[ |33p (BGR), [|25| (RBC N F = 
2 + 1), [|19|^] (RBC N F = 2), [Q (LHPC), |$5j (ETMC) and [[24]] (QCDSF). The experimental result is the 
PDG 2008 value fR8|). The discrepancy between the lattice calculations and the experimental measurement 



and the scatter among the lattice calculations are both smaller than for shown in Fig. ^[ 



Fig. ^ by examining the lighter pion mass calculations of both QCDSF and ETMC. It is possible 
that the volumes for these calculations at m K = 310 MeV or 260 MeV may still be too small to 
see the asymptotic volume dependence. We can also consult chiral perturbation theory, which does 
allow for this sort of behavior; however, given that gA is essentially flat for the largest volumes, one 
might reasonably question the use of chiral perturbation theory at these pion masses. 

In order to attempt to estimate the infinite volume limits for the lattice calculations presented 
here, I simply require m n L > 6. This may turn out to be an excessive requirement for low m K , 
which would actually be good news regarding finite size effects, but it is clearly required for heavier 



pion masses. The lattice results satisfying this are shown in Fig. 10. As is clear from the figure, 
this restriction is very severe, however, the resulting lattice calculations show a strong level of 
agreement amongst themselves. The results are still lower than the experimental measurement, but, 
unlike for (x) u _^, only a mild curvature is required to reconcile the current calculations with the 
physical limit. In particular, notice that, with one exception, there is no systematic shift between the 
various collaborations. As mentioned earlier, the renormalization of gA is generically easier than 
that of (x) u -d> an d this lends a bit more support to the hypothesis that differences in renormalization 
may be driving part of the variation of (x)„_d in Fig. ^. The one small exception is the result of 
QCDSF which is just slightly low compared to all the other calculations. However, this may be 
consistent with a small discrepancy in f n , which renormalizes with the same factor as gA, that is 
also present for QCDSF | 
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Figure 9: Finite size effects in g^. The m n L dependence is shown for several calculations. The results at the 
heaviest pion masses indicate a significant finite size effect in g^, suggesting that m K L of 6 may be necessary. 
But the light m K results of QCDSF and ETMC may in fact be showing a weakening of the finite size effects 
in gA as measured in terms of m K L. The lattice results and experimental measurement are the same as in 
Fig. H and the curves are simple scaling fits to guide the eye. 



4. Nucleon Charge 



The charge of the nucleon, or any hadron, again appears to be an essentially trivial topic. 
There is no doubt that the proton has one unit of charge and that the neutron is, well, neutral. 
The certainty in our understanding of the nucleon's charge arises from the flavor symmetries of 
QCD, however, the distribution of this charge in the spatial degrees of freedom of the nucleon's 
constituents is not prescribed by symmetry and provides another probe of the dynamics of QCD 
within the nucleon. This idea leads to a broad range of calculations on the lattice and quickly 
involves the generalized partem distributions, but to avoid excessive complications and to explain 



the concepts in the simplest setting, I will discuss just the neutron charge [ |37| ] 



4.1 Neutron Transverse Charge Distribution 

Much of the interest in form factors originates with the interpretations that we can assign to 
them. There is now a rich field theoretic discussion dedicated to this issue, but in order to maintain 
our footing and retain a strong connection to the experimental measurements, it is valuable to 
remember that ultimately, the form factors parametrize the QCD contributions to cross sections, 
independent of any interpretation we attach to them. 

^=(^)M„, {F ' + ^ + 2T(F|+F2)2 ' a " 2(e/2,} 
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Figure 10: Results for g& with m n L > 6. The lattice results from Fig. |]that satisfy m K L > 6 are shown. All 
the lattice calculations, but one, are in statistical agreement. The results still underestimate the experimental 
value, but the discrepancy is much less severe than in (x)„_ t /. The QCDSF results are just slightly low. 
This may indicate some small systematic error in the renormalization or it may simply illustrate the lattice 
artifacts that are present in all the calculations. 



Here F\ (q 2 ) and /^(g 2 ) are the nucleon form factors and z = q 2 / (4m 2 ). For some in our commu- 
nity, this is where the discussion ends. That is, of course, a perfectly valid point of view. However, 
the understanding provided to us by the generalized parton distributions and the relationship to 
transverse quark distributions [38] allows us to go further. 

Historically the form factors were interpreted as charge and magnetization densities. In par- 
ticular the charge density was constructed from Ge = F\ — 1F2 as 

J (2ny ^J m 2 + ji 

This is plotted for the neutron in Fig. [TT| There we see that the charge distribution is negative for 
large r. This is usually understood in terms of the pion cloud generated by n — > p%~ fluctuations. 
Overall neutrality of the neutron forces the core to be positive leading to a negative charge radius. 

But this picture is spoiled by relativistic corrections that enter as (m(r))~ l , where (r) is some 
relevant length scale for the charge distribution. These corrections would vanish in the non- 
relativistic limit m — > oo but are substantial for the nucleon. The introduction of the transverse 
distributions in the infinite-momentum frame (IMF) eliminates these corrections, which enter now 



as (p z (r)) 1 with p z — > oo defining the IMF. The transverse distributions [38] are given by 

This is shown for the neutron in Fig. O. The distribution still has a negative tail for large r, but 
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Figure 11: Naive 3d neutron charge distribution. The 3d Fourier transform of the Ge form factor gives the 
rest frame charge distribution of the nucleon up to relativistic corrections. The negative distribution at large 
r is interpreted as the pion cloud arising from the n — > pn~ fluctuations. This gives rise to a negative charge 
radius for the neutron. The curves were generated by numerically integrating the Kelly 2004 parametrization 



now we find a negative core. Furthermore, the charge radius is positive for the transverse charge 
distribution. Understanding this new view of the neutron is a challenge onto itself. Additionally, 
the introduction of the transverse coordinates in a field-theoretically clean manner opens many new 
avenues for investigation of hadron structure from lattice QCD. 

4.2 Lattice Calculation of (r 2 ) u d 

As for the nucleon momentum and spin structure, I again introduce a benchmark observable 
that allows us to summarize the status of lattice calculations of form factors. Here I choose the 
isovector charge radius, (r 2 )„_^. As before, the flavor structure is u — d, so disconnected diagrams 
again vanish. Furthermore, this observable enters matrix elements of the vector current. These two 
facts combine to make (r 2 ) H _^ an accurate quantity to calculate on the lattice. Additionally, all the 
lattice calculations can now use an exactly conserved vector current. This eliminates the need for 
any renormalization and simplifies the calculation significantly. 

The form of the nucleon's electromagnetic interaction is given by 



of the nucleon form factors 



p,s)=u(p>,s>){ 7 »F l (q 2 ) + 




2m 



■F 2 (q 2 )}u(p,s) 



and the charge radius is then given by 



dF, 



u—d 




6 



dq 2 
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Figure 12: Transverse 2d neutron charge distribution. The 2d Fourier transform of the F\ form factor results 
in the transverse charge distribution of the nucleon. Relativistic corrections are absent but the distribution 
must be understood in the infinite momentum frame. Notice the negative core, in contrast to Fig. [TT| Also 
the charge radius is now positive in this frame. The distribution was constructed from the same form factor 
parametrization in Fig. [HJ 



The factor of 6 is the conventional choice, but a factor of 4 would instead give the correct transverse 
charge radius in the IMF. To avoid confusion, we stick to the conventional choice, but we should 
keep in mind the theoretically clean interpretation that is available for the transverse distribution. 

In Fig. 13 I gather all the full QCD results for (r 2 ) u -d- It is immediately evident that all 
the lattice calculations show a remarkable agreement for this observable. The strong scatter that 
is present for (x} u -d and to a lesser extent gA is clearly absent for (r 2 )„_d. This is particularly 
surprising given the fact that this quantity is expected to diverge in the chiral limit as indicated by 
the curve in Fig. 13. We might even expect large finite size effects in such a quantity. However, 
we must keep in mind that {r 2 ) u -d is calculated from the slope of Fi(q 2 ) at q 2 = and that the 
momentum quantization imposed by a finite volume forces us to extract this slope using the first 
non-zero value of q 2 , which satisfies qL « 271. We may, in fact, be extrapolating from a region 
with weaker L dependence through the region with stronger L dependence. Ultimately, whatever 
the explanation, we find a near absence of volume dependence for (r 2 )„_^. 

Lacking any real guidance on the necessary volumes for (r 2 )„_</, I simply follow the conven- 
tional wisdom and take m K L > 4. These lattice results are shown in Fig. 14 . An optimist may find 
a mild increase, but this would be hard to defend strongly. Additionally, this observable requires 
much more than a mild increase to agree with the experimental value at the physical point. But it 
does appear that all the calculations seem to agree. Given that {r 2 ) u -d requires no renormalization, 
this is consistent with the speculation suggesting that renormalization may be partially culpable for 
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Figure 13: World's dynamical results for (r 2 )^. The lattice calculations are Q (RBC Np = 2 + 1), 
[[9[ |§ (RBC N F = 2), [0 (LHPC/MILC), [|IJ] (LHPC/RBC), [||] (ETMC) and @] (QCDSF). The 
experimental result comes from the Kelly 2004 parametrization of the form factors [pSfl. The leading order 
heavy baryon chiral perturbation theory prediction (LO HB^PT) contains no unknown low energy constants 
other than (r 2 )]/^ and is logarithmically divergent as m K — » 0. One might, therefore, expect substantial finite 
size effects in this quantity, but this does not appear to be the case. 

the systematic variation in {x) u -d discussed earlier. Despite the fact that this observable still shows 
a significant discrepancy with the experimental result, the agreement amongst all the lattice groups 
is encouraging. 

5. Conclusions and Outlook 

In order to illustrate the status and promise of lattice calculations of hadron structure, I chose 
to review three benchmark observables that represent the broad range of hadronic physics that we 
can address on the lattice: the momentum fraction, axial charge and charge radius of the nucleon. 
Nearly all the lattice groups have pushed the pion mass towards 300 MeV or lower and some are 
now breaking through 200 MeV. This is an enormous accomplishment, but unfortunately most 
hadronic observables are failing to show the long-sought chiral curvature. Independent of the 
expectations of chiral perturbation theory, the three benchmark observables still do not show a 
convincing approach to the experimental results. 

The obvious conclusion is that yet lighter, and maybe even physical, pion masses will be 
required. This may be the case, but we should not do this at the expense of other systematic errors. 
Finite size effects must remain a concern as we lower m n . Both (x) u -d and gA demonstrate that the 
volume dependence can have a non-trivial dependence on m n and our conventional rule-of-thumb 
that m n L « 4 may not necessarily be sufficient for all quantities. 
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Figure 14: Large volume results for {r 2 )){_ d - Due to momentum quantization at finite volume, one can 
not easily calculate the form factors at the low Q 2 needed to determine {r 2 )\j^ d while also working at small 
volume. Lacking any detailed knowledge of the finite size effects in (r 2 ) l J 2 d , only the results from Fig. [l3] 
with m n L > 4 are shown. There is general agreement among all the lattice groups and a mild increase as 
m K is lowered. However (r 2 )]/ 2 d is not yet rising nearly as fast as it needs to match onto the experimental 
measurement. 



Lattice artifacts are still a concern, if only because so few groups have seriously examined 
them. The agreement of the groups for (r^) u -d might indicate weak artifacts for some quantities, 
but the systematic shifts between the groups for (x) u -d is a warning that renormalization, an essen- 
tial part of the continuum limit for most quantities of interest, may be partially responsible for these 
problems. The role of renormalization is further implicated by the general agreement of the groups 
for the renormalization free ratio of {x) u -d/ '{ x Yu-d- m ^ act ' tne progression from {x) u -d (requir- 
ing a scale-dependent renormalization) to gA (requiring only a finite renormalization) to (r 2 ) H _^ 
(requiring no renormalization) will be a strong test of our control of renormalization issues. 

These systematic effects will be a part of any critical review of lattice calculations, but we 
should remember the rich QCD physics that provides the motivation for these calculations. Each of 
the benchmark observables is a stand-in for a broad range of experimentally relevant and intellec- 
tually interesting observables: the decomposition of momentum among the nucleon's constituents, 
the nature of the nucleon's spin and the distribution of the nucleon's charge. Apart from these 
observables, Lattice 2009 saw a sizeable effort to study strangeness in the nucleon, results for the 
pion, delta and other hadrons, and a range of exploratory calculations. Clean control of systematics 
and a continued push toward the physical pion mass will ultimately allow detailed calculations that 
will complement the experimental programs and advance our insight into hadron structure. 
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